Connected-Sea Partons 



Keh-Fei Liu, 1 Wen-Chen Chang, 2 Hai-Yang Cheng, 2 and Jen-Chieh Peng 3 

1 Department of Physics and Astronomy, University of Kentucky, Lexington, Kentucky 40506, USA 

2 Institute of Physics, Academia Sinica, Taipei 11529, Taiwan 
3 Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA 

According to the path-integral formalism of the hadronic tensor, the nucleon sea contains two distinct compo- 
nents called connected sea (CS) and disconnected sea (DS). We discuss how the CS and DS are accessed in the 
lattice QCD calculation of the moments of the parton distributions. We show that the CS and DS components 
for u(x) + d(x) can be extracted by using recent data on the strangeness parton distribution, the CT10 global 
fit, and the lattice result of the ratio of the strange to u(d) moments in the disconnected insertion. The extracted 
CS and DS for u(x) + d(x) have distinct Bjorken x dependence in qualitative agreement with expectation. The 
analysis also shows that the momentum fraction of the u(x) + d(x) is about equally divided between CS and 
DS at Q 2 — 2.5GeV 2 . Implications on future global analysis for parton distributions are presented. 
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There have been a number of developments in the under- 
standing of the flavor content of the nucleon sea, such as 
the observation of the light-quark sea difference between d 
and u in Deep Inelastic Scattering (DIS) flj and Drell-Yan 
processes 0, the extraction of strange quark content s + s 
from semi-inclusive DIS [3], and the lattice QCD calcula- 
tions of sea quark contributions to nucleon orbital angular mo- 
menta [ 4 1 . Evidence for the existence of intrinsic sea Q of the 
light quarks has also been reported [6 1. 

Many theoretical models, including the meson cloud 
model, have been suggested for describing the flavor structure 
of the nucleon sea |7 1. In order to gain new insights on the ori- 
gins of the flavor content of the nucleon sea, it is important to 
note that, according to the path-integral formalism, there are 
two distinct sources for nucleon sea, namely, the connected 
sea (CS) and the disconnected sea (DS) MM- The CS and 
DS are expected to have different shapes in their Bjorken-x 
distribution, as well as distinct quark-flavor dependence. The 
first direct experimental evidence for the existence of CS came 
from the observation of large difference in the u(x) and d(x) 
distributions [8]. In this paper we show that the two distinct 
contributions (CS and DS) to the u(x) + d(x) can be sepa- 
rated based on existing experimental data and input from lat- 
tice QCD calculation. 

The existence of the connected sea and disconnected 
sea can be illustrated in the path-integral formalism of the 
hadronic tensor. In the Euclidean path-integral formalism of 
the hadronic tensor W^, there are three gauge invariant and 
topologically distinct diagrams, as shown in Fig. 1. The var- 
ious lines in Fig. 1 represent the quark propagators from the 
source of the nucleon interpolation field at time t = to the 
sink time at t and the currents are inserted at t\ and t^. 



We first note that Fig. 1(b) where the quarks propagate 
backward in time between t\ and t%, corresponds to the 
connected-sea anti-partons u cs and d cs . In contrast, the for- 



ward propagating quarks in Fig. 1(a) correspond to valence 
and CS partons u v+cs and d v+cs , where valence is defined as 

q v = q v+cs _ -cs and q cs^ 

gives the DS q ds and q ds for q 



both forward and backward propagating quarks. The nomen- 
clature of connected and disconnected seas follows those in 
the time-ordered perturbation theory - CS is the higher Fock- 
state component in the Z-graph where the quark lines are con- 
nected and the DS corresponds to the vacuum polarization. 
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FIG. 1 : Three gauge invariant and topologically distinct diagrams in 
the Euclidean path-integral formalism of the nucleon hadronic tensor 
in the large momentum frame. In between the currents at ti and ti, 
the parton degrees of freedom are (a) the valence and CS partons 
q v+ca , (b) the CS anti-partons q ca , and (c) the DS partons q da and 
anti-partons q da with q = it, d, s, and c. Only u and d are present in 
(a) and (b). 

It is clear from Fig. 1 that the two sources of the sea quarks, 
CS and DS, have interesting quark-flavor dependence. For ex- 
ample, while u and d have both CS and DS, s and c have only 
DS. The small mass difference between the u and d quarks 
implies that the DS cannot account for the large d/u differ- 
ence observed in the DIS and Drell-Yan experiments. Rather, 
this difference must originate primarily from the CS diagram 
of Fig. 1(b) due to the fact that there are two u-valence quarks 
but only one d. The absence of the CS component for the 
strange and charm quarks also implies that any difference be- 
tween s(x) and s(x) (or c(x) and c(x)) distributions, as pre- 
dicted in meson cloud 1 10] and intrinsic sea ifTTI models, must 
come from the DS diagram of Fig. 1 (c). The classification of 
parton distribution functions in terms of flavor, CS and DS is 
given in Table [I] 

The CS and DS are also expected to have distinct distribu- 
tions at the small-x region. Since there is only reggeon ex- 
change for the flavor non-singlet valence and CS, the valence 
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FIG. 2: The three-point functions after the short-distance expansion 
of the hadronic tensor from Fig. 1. (a) The connected insertion 
(CI) is derived from Fig. |l(a)| and Fig. |l(b)| (b) The disconnected 
insertion (DI) originates from Fig. 1(c) Og are local operators which 



are the same as derived from OPE 



TABLE I: Classification of PDF in the nucleon for different flavors. 



Valence and Connected Sea 


u u+ca (x) 


u cs {x) 


d v+cs (x) 


d cs (x) 


Disconnected Sea 


u da (x) + u d3 (x) 


d ds {x) + d d3 (x) 


s(x) + s(x) 


c(x) + c(x) 



and CS partons is q v+cs (x), q cs (x) — > q oc x 1//2 at small x. 
For the DS partons, there is flavor-singlet pomeron exchange, 
thus its small-x behavior is q ds {x) 1 q ds (x) — > q oc x _1 . In 
addition, there are also meson cloud contributions which are 
prominent in the medium- a; range. Partons with these differ- 
ent small-x behaviors are considered as extrinsic and intrin- 
sic distributions for charm [5| and the light quarks J6). The 
distinct flavor and x-dependence of CS and DS remain to be 
checked experimentally. 

Under the short-distance expansion of the hadronic tensor 
between the current insertions in the path-integral formalism 
(N.B. This corresponds to operator product expansion (OPE) 
in the canonical formalism), Fig. 1(a) and Fig. 1(b) become 



the connected insertions (CI) in Fig. 2(a) for a series of local 
operators y] n O™ in the three-point functions from which the 
nucleon matrix elements for the moments of the CI are ob- 
tained. Here the flavor q = u,d are the valence flavors from 
the interpolation field. By the same token, the disconnected 



four-point functions in Fig. 1(c) become the disconnected in- 
sertions (DI) in Fig. |2(b)| for the three-point functions to ob- 
tain the DI moments. Here q = u,d,s,c are the DS flavors 
in the DI. The main advantage of the path-integral formalism 
over the canonical formalism is that the parton degrees of free- 
dom are tied to the topology of the quark skeleton diagrams 
in Figs. [T(a)l [Kb)] and [T(c)1 so that the CS and the DS can be 
separated. Lattice QCD can access these three-point functions 
for the CI and DI which separately contain the CS and DS and 
calculations of the moments of the unpolarized and polarized 
PDFs for the quarks [12] and glue [4| have been carried out. 



Unfortunately, lattice calculations cannot calculate the parton 
x-distributions directly [9|, only moments are accessible. 

To delineate the flavor and x dependence of the parton dis- 
tributions, it is important to have the CS and DS separated in 
the global fitting as they evolve differently in Q 2 |9|. Once 
they are separated in one Q 2 , they will remain separated so 
that they can be used to fit experiments or make predictions 
at other Q 2 . While the difference of u cs (x) and d cs (x) is 
obtained from the E866 Drell-Yan and HERMES semi- 
inclusive lfL3l measurements of d(x) — u(x), there is not yet 
an established way to directly obtain u cs (x) + d cs (x) from 
experiments. We shall show how to achieve this with a com- 
bination of experimental results and a lattice calculation of (x) 
for the DI in Fig. [2(b)] 
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FIG. 3: The strange quark PDF from HERMES semi-inclusive DIS 
experiment of kaon production on deuteron. It is compared to CT10 
results. 

Recent HERMES semi-inclusive DIS experiment of kaon 
production on deuteron has produced the strangeness par- 
ton distribution function s(x)+s(x) at Q 2 = 2.5 GeV 2 which 
is shown in Fig. [3] One notable feature is that the values at 
medium x around 0. 1 are quite different from those obtained 
from the global fit of CT10 [ 14- 1 (also drawn in Fig. [3J, which 
did not include the HERMES data in the fit. Since s and s 
are entirely due to DS, the HERMES s(x) + s(x) data pro- 
vide valuable information on the shape of the x-distribution 
for DS, which is not available from the lattice calculation. 

We can now proceed to separate the CS and DS compo- 
nents of the u + d sea with the following approach. First, 
we shall make the plausible ansatz that the distribution of 
u ds (x) + d ds (x) is proportional to that of s(x) + s(x) and the 
proportionality is j^, i.e., u ds (x) + d ds (x) — j^(s(x) + s(x)). 
As discussed below, a recent lattice calculation has obtained 
R = 0.857 ± 0.040. We can then extract u cs (x) + d cs (x) 
from the difference of the CT10 result on u(x) + d(x) and the 
HERMES data weighted with ^, 

u cs {x) + d cs (x) = u(x) + d(x) - 4O0*0 + s(x)) (1) 

it 
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at Q 2 = 2.5 GeV 2 . The u(x) + d(x) in Eq. 1 can be taken 
from the recent CT10 PDF. 

A recent lattice calculation lfT31 of the momentum fraction 
(x) for the strange and u(d) in the DI was carried out with 
2 + 1 -flavor dynamical fermion gauge configurations with the 
improved gauge and Wilson fermion (clover) actions. The 
ratio of the momentum fraction between the strange and the 
u(d) in the DI, where multiplicative renormalization constants 
and some systematic errors cancel, is extrapolated to the chi- 
ral limit with the lowest pion mass at 600 MeV and the result 



R 



(X) 



s + s 



(x) u+a (DI) (x) udB+dds 



0.857(40). (2) 



The second equality is based on the premises of isospin sym- 
metry of the DS, i.e., u ds — d ds and that the parton-antiparton 
difference is negligible in the DS, i.e., u ds — u ds . 

It is interesting to note that the ratio R in Eq. d2l) is close 
to unity, much larger than the globally fitted ratio of ) / J+g _ ~ 
0.5 as evidenced in the CT10 results in Fig. [3] The difference 
is due to the fact that u + d has an additional contribution 
from CS which the strange partons do not have. Moreover, the 
ratio R = 0.857 is very close to that of the CT10 fit at small 
x (e.g., x < 2 x 10~ 2 ). This is consistent with the recent 
ATLAS measurement of the inclusive W and Z productions, 
where the strange-to-down quark ratio was determined to be 
1.00±g;|| at x = 0.023 and Q 2 = 1.9 GeV 2 ED. Since 
partons at small x are dominated by the DS, this shows that 
both the ratio of (s(x) + s(x))/(u ds (x) + d ds {x)) of CT10 
(and ATLAS) at small x and the ratio of their second moments 
from the lattice are practically the same. This supports the 
ansatz that u ds (x) + d ds (x) = j^(s(x) + s(x)). 

We plot the distribution function evaluated with Eq. ([T}, 
multiplied by the momentum fraction, i.e., x(u(x) + d(x) — 
j^(s(x) + s(x)) in Fig. [4] together with x(d(x) — u(x)) from 
E866 Drell-Yan measurement J2] at Q 2 = 54 GeV 2 and from 
SIDIS HERMES measurement JB) at (Q 2 ) = 2.3 GeV 2 . 
We see that x(u cs (x) + d cs (x)) from Eq. ([IJ is peaked at 
medium x ~ 0.1, the same way as x(d(x) — u(x)) from 
E866 and HERMES. This is consistent with the expectation 
that the small-a; of CS, like the valence, behaves as x~ x l 2 as 
we alluded to earlier; so that, when CS is multiplied with x, 
it would be peaked at medium x, in contrast to that of the 
DS, e.g., x(s(x) + s(s)) in Fig. [3] Furthermore, we note that 
x(u cs (x) + d cs (x)) is generally larger than x(d(x) — u(x)) 
in this re-range as it should and is larger by a factor ~ 4 at the 
peak. 

We also plot x(u(x) + d(x) — j^(s(x) + s(x)), x{u ds {x) + 
d ds (x) = jtx(s(x) + s(x)) and x(u(x) + d(x)) from CT10 
in Fig. [5] to show that the CS and DS have very different re- 
dependence. The different shapes of CS and DS are in good 
agreement with the expectation discussed earlier. This agree- 
ment lends support for the approach we adopted. It is inter- 
esting to note that should a very different value of R be used, 
the x-dependence of CS and DS would no longer agree with 
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FIG. 4: x(d ca (x) + u ca (x)) obtained from Eq. |T} is plotted together 
with x(d(x) — u(x)) from E866 Drell-Yan experiment [j2] and from 
SIDIS HERMES experiment (TJ. 



expectation. In particular, if R were appreciably larger than 
the present value of 0.857, Eq. (fill would lead to a CS whose 
small x behavior would be ~ x which is inconsistent with 
the fact that CS is from the connected insertion. In this case, 
x(u cs (x) + d cs (x)) would approach a constant as x —> as 
opposed to zero, as illustrated in Fig. [5] On the other hand, if 
R were appreciably smaller than 0.857, the CS from Eq. ([T]) 
would turn out to be negative and this would not be commen- 
surate with the probability interpretation for CS. 
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FIG. 5: x(u ca (x) + d ca (x)) obtained from Eq. |T} is plotted together 
with x(u(x)+d(x)) from CT10 and j^x(s(x) + s(x)) which is taken 
tobex(u da {x) + d ds {x)). 

We can also calculate the momentum fractions carried by 
the CS and DS of u + d as follows. First, the total con- 
tributions of u + d, i.e., (a;)j +a = L dxx(d(x) + u(x)) for 
the global fittings of CT10, CTEQ6 Q7) and MSTW08 HI 
are listed in Table II. Here, the uncertainty of a specific PDF 
function is estimated by varying the range of the calculated 
quantity with input from the alternative eigenvectors of the 
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TABLE II: Values of various moments using the HERMES data, the 
lattice QCD result, and three different PDFs at Q 2 = 2.5 GeV 2 . 





CT10 


CTEQ6 


MSTW08 


{ x )d+u 


0.0639(14) 


0.0614(14) 


0.0690(11) 


( X ) d CB +u cs 


0.0294(54) 


0.0281(54) 


0.0347(53) 


i X ) d dB +u ds 


0.0344(52) 


0.0332(52) 


0.0342(52) 


(*) d ds +a ds 


0.86(29) 


0.85(30) 


1.02(32) 



PDF and is quoted in the parentheses. The CS contributions 
(x) +il c S in Table|lT|are then obtained from Eq.[T|with inputs 
of d(x)+u(x) from the different PDFs and integrated over the 
range 0.025 < x < 0.48 plus a small contribution outside this 
range which is estimated to be 2% from the fractional contri- 
bution to x(d(x) — u(x)) outside this range based on CT10. 
The total uncertainty includes the contributions from HER- 
MES data, the lattice calculation of R, the PDF, and the un- 
measured x region. The corresponding (x) gds +U da in Table|n] 
is taken to be the difference between the respective {x)^ +n 

(x) -cs -cs 

and (x) jesses. Finally, we give the ratio ) ! d _ +c . It is in- 

- \ x ) d ds + ^ds 

teresting that this ratio is close to unity, showing that the mo- 
mentum fraction of d + u is about equally divided between 
the CS and the DS at this low Q 2 . Future lattice calcula- 
tions could provide a direct check of the second moment of 
u ds (x) + d ds (x) shown in Table II. 

In order to gain deeper and more precise understanding of 
the PDF in terms of their flavor, x and Q 2 dependence, it is 
essential to have the CS and DS separately accommodated in 
the extended evolution equations. Only then will they remain 
separated at different Q 2 to facilitate global fitting, as they 
evolve differently with the CS evolving like the valence J5). 
This will have an impact on the gluon distribution as well. 
Furthermore, only with CS and DS separated, will one be 
able to address the flavor dependence, i.e., u ds ^ d da and 
parton-antiparton difference of the u and d partons in the DS 
and check the validity of the ansatz that u (x) + d ds (x) is 
proportional to s{x) + s(x). As the lattice calculations are 
getting more refined when the physical pion mass, continuum 
limit, and large volume limit are approached, they could serve 
as valuable constraints for the parton moments. In particu- 
lar, higher moments in CI will help separate the valence and 
CS parton distributions in the global analysis and the DI cal- 
culation of the fourth moment for the strange and u/d (i.e., 
(x 3 ) u (DI) = J dxx 3 (u ds (x) + u ds (x))) can be used to 
gauge how good the proportionality assumption is about their 
distributions. 

In summary, we have shown that there are two sources for 
the sea partons, the CS and DS, based on the path-integral for- 
malism of the hadronic tensor. While the u and d have both 
CS and DS contributions, the s and c partons are from the 



DS only. We also expect that the CS and DS have different 
x distributions. These different flavor and x dependence offer 
the possibility of disentangling the CS from the DS. We first 
show that the expectation for the dominance of DS at small 
x is supported by the good agreement between the lattice cal- 
culation of R = = 0.857(40) and the ratio of 
x(s(x) + s(x)) from the HERMES data to x(u{x) + d(x)) 
from CT10 at small x (e.g., x < 2 x 10~ 2 ) as shown in Fig|3] 
Given this agreement, we show how the HERMES data on 
strangeness parton distributions, the lattice calculation of R, 
and the CT10 global fit of d + u, can determine the separate 
CS and DS contributions in u + d. We stress that the ansatz of 
the proportionality between u ds (x) + d ds {x) and s(x) + s(x) 
should be checked with lattice calculation of R with light dy- 
namical fermions as well as the ratio for (x 3 ) in the DI. Future 
global analysis of PDF should have CS and DS separated in 
the fitting and in the evolution equations. 
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